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The exact expression derived by Bougourzi, Couture, and 
Kacir for the 2-spinon contribution to the dynamic spin struc- 
ture factor S zz (q,u)) of the one-dimensional S=l/2 Heisen- 
berg antiferromagnet at T = is evaluated for direct com- 
parison with finite-chain transition rates (TV < 28) and an 
approximate analytical result previously inferred from finite- 
TV data, sum rules, and Bethe-ansatz calculations. The 2- 
spinon excitations account for 72.89% of the total intensity 
in S zz (q,uj). The singularity structure of the exact result 
is determined analytically and its spectral-weight distribu- 
tion evaluated numerically over the entire range of the 2- 
spinon continuum. The leading singularities of the frequency- 
dependent spin autocorrelation function, static spin structure 
factor, and g-dependent susceptibility are determined via sum 
rules. 



I. INTRODUCTION 

Notwithstanding the fact that BetheS found the key 
that solves the one-dimensional (ID) 5=1/2 Heisenberg 
model, 



ff = jJ]S|-S H 



(1.1) 



i=i 



as early as 1931, the emergence of explicit results for 
various physical-Quantities (ground-state energy,0 excL 
tation spectrum J magnetization curve, susceptibilityjj'tl 
thermodynamic^]) was slow at first and then faster since 
around 1960. Interest in this model began to spread far 
and wide when the first compounds with quasi-lD mag- 
netic properties were synthesized and investigated exper- 
imentally. 

However, the dynamics of th e ID Heisenberg antiferro- 
magnet, i.e. Hamiltonian ( |l.l| ) with J > 0, has remained 
elusive to any rigorous approach during all those years. 
An exact result for the dynamic spin structure factor 



S zz (q,w) 



l,n 



dte luJt {St(t)S : 



l+nl 



(1.2) 



in particular, would have been of great value for the in- 
terpretation of a host of experimental data.Q 

Significant progress in the understanding of the T = 
dynamics resulted from the observation^ that almost all 



the spectral weight in S zz (q,u>) is carried by a special 
class of Bethe-ansatz solutions with excitation energies 
(in units of J henceforth) 



/ \ . q (q q m 

uj m (q) =7rsm-cos - — 



(1.3) 



for N — > oo, 0<g<7r, 0< q m < q. They form a two- 
parameter continuum in the (g, w)-plane bounded by the 
branches 

UL(q) = | sing, uiu(q) =7rsin-. (1.4) 

These excitations were later named 2-spinon states. 
Thjeir density of states (after rescaling by a factor 2tt/N) 



D{q,L0) 



LO L (q))Q(u u (q) - J) 

V^uiq) - 



(1.5) 



The T = dynamic spin structure factor for a finite 
system with even TV and periodic boundary conditions 
can be written in the form 



S zz (q,uj) = 2n} j M\5(u - 



(1.6) 



where M x = \ (G\S z q \X)\ 2 with S z q = N~ x / 2 e^5f are 
the transition rates between the singlet (St = 0) ground 
state \G) and the triplet (St = 1) states |A) with finite- 
TV excitation energies u>\. Among them are the N(N + 
2)/8 2-spinon excitations, which contribute most of the 
spectral weight. 

The finite-chain analysis of Ref. || suggested that the 
scaled transition rates NM\ vary smoothly with q and uj. 
The consequence could be that the exact 2-spinon part of 
the dynamic structure factor is expressible, in the limit 
TV — > oo, as a product 



Si 2 J(q,Lu) = M(q,L0)D(q,u), 



(1.7) 



with a smooth transition-rate function M(q,uj), toward 
which the scaled finite- TV transition rates converge. Tip*] 
scenario is indeed realized in the related XX modelJat£l 
where the 2-spinon density of states is given by ( |l.5| ) 
with modified spectral boundaries, and the transition- 
rate function is a constant .B 
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II. APPROXIMATE TRANSITION RATES 



In the Heisenberg model (1.1), the finite- N data for the 
2-spinon matrix elements indicate that M(q, uj) diverges 
at u> = and vanishes at u) = uju(q). In Ref. || the 

expression 



MW(?,( 



(2.1) 



for the 2-spinon transition-rate function was proposed 
on the basis of this observation and the following three 
requirements: The resulting (approximate) 2-spinon dy- 
namic structure factor S z t\q,uj) = M (a \q,uj)D(q,uj) 
must produce, for q — tt, the correct infrared 
exponent O'EHl Furthermore, it must produce the copjject 
g-dependence of the known first frequency moment jp t3 



Ki(q) 



2^ 



wS zz (q, uj) 



2Eg 
' 3N 



(1-cosg), (2.2) 



where Eq — — iV(ln2 — 1/4) is the ground-state energyE 
and, via the sum rule 



1 



du> 



X (q) = - / —S„{g,u), (2.3) 
n J uj 
o 

the correct value for the direct susceptibilitvlp] x(0) = 
I/71- 2 . The resulting approximate expression JiJ 



Qjojj- LO L (q))Q(uj u (q) - u)) 
y/u 2 -u)l{q) 



(2.4) 



for the 2-spinon dynamic structure factor has been used 
quite frequently for the interpretation of inelastic neutron 
scattering measurements on a number of quasi- ID anti- 
ferromagnets at low temperature,Q and for comparisons 
with the results of various computational studies .EjtZl 

It is interesting to note in this context that the ex- 
act dynamic structure factor S zz (q,u>) of the Haldapc- 
Shastry model has a structure very similar to (2.4) Jia In 
that model, as in the XX model, all the spectral weight 
of S zz (q,ui) is carried by the 2-spinon excitations. 



III. EXACT TRANSITION RATES 



A detailed assessment of the merits and limitations of 



the result ( j2.4| ) has become possible only recently through 
a remarkable new development. By approaches based on 
the concept of infinite-dimensional symmetries which had 
been developed in the context of siring theory, conformal 
field theory and quantum groupstj Bougourzi, Couture, 
and KacircJ were able to derive the exact expression for 
the 2-spinon transition-rate function in the formed 



M( q> u) = ~e- I V 



where t = 2(/3i — /?2)/7r and 
I(t) = 



, cosh(2x) cos(a;t) — 1 „ 
dx , , — , — ; e 



xsinh(2a;) coshx 



2 cosh (3i 2 cosh (3% ' 
q= — cot (sinh/3i) — cot (sinh/^). 



(3.1) 

(3.2) 

(3.3a) 
(3.3b) 



By solving Eqs. (gj) we can express the auxiliary vari- 
able t as a function of the two physical variables q, uj: 



4 j ujl(q)-ujl(q) 
UJ 2 -uj 2 L (q) 



t = —cosh 

7T 



(3.4) 



For the numerical evaluation of (3.1 ) we separate the sin- 
gular part from the integral (3.2): 



/(*) = - Jo - In ^isinh 2 + h(t), 

where 

h(t) = Ci(t) + fx(t) - f 2 (t), 

oo 1 

j. . . f dx cos(xt) j . . f dx cos(xt) 
J x cosh 2 x ' J x coth 2 x ' 

1 

Jo = 7 + /i(0)-/a(0)= 0.3677103.... 
A series expansion of s(t) = [lni + C — h(t)]/2, 

oo 

sin 2 (xi/2) 



(3.5) 

(3.6a) 
(3.6b) 
(3.6c) 



s(t) = / dx 2 — 

' x cosh x 



= ^(-l) ro mm 1 



m—1 



4m z 



C = e Io /2 = 0.72221 



(3.7) 



brings (|3.1[) with t from (3.4) into closed form 



M{q,uj) = Ctsinh 



n 



{l + [t/(4m-2)] 2 } 
{1 + [t/4m] 2 } 2m 



2i 2m— 1 



(3.8) 



Th e exact 2-spinon part of S z z (q, uj), i.e. the function 
(1.7) with the den sity of states ( [L5| ) and the t rans ition- 
rate function (3^) evaluated numerically via (3J3) with 
( |3.4[ ) is plotte d in Fig. 1(a). For comparison, the approx- 
imate result (2.4) is shown in Fig. 1(b). The two results 
look very similar, yet there are subtle differences, which 
may not matter for most experimental comparisons but 
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are important for comparisons with other theoretical re- 
sults. 

Both expressions diverge at the lower spectral bound- 
ary u>l(q)- At the upper boundary uju{q), S z a z \q,uj) has 
discontinuity, whereas S zz {q, oj) approaches zero contin- 
uously over a rounded shoulder. 

The structure of the exact transition rate function (3.1) 
lends itself naturally to be factorized into the approx- 
imate function ( |2.l| ) and a correction which accounts 
for the modified singularities at the boundaries of the 
2-spinon continuum: 



M(q,u) = M { - a \q,oj)^/Ct/2e h ^' 2 



(3.9) 



IV. SINGULARITIES AT u L (q) AND oj v {q) 

What is the precise nature of the leading singular- 
ity in the transition-rate function M(q,uj) and in the 

2-spinon dynamic structure factor S zz (q,uj) at the spec- 
tral boundaries uJu(q) and wl(q), and how do these sin- 
gularities compare with those of the approximate results 
S^fewjand M^J^u;)? The answer is obtained by in- 
serting (3.5) into (3.1), evaluating the leading term for 
t — > and t — » co, respectively, and inserting (3.4) ex- 
panded accordingly. 

At tou(q) the transition-rate function is thus found to 
approach zero linearly, 



M(q,u) 



8C 



uu{q) 



TT ujfj(q) - LJ 2 L (q) 



[uu(q) - u], (4.1) 



which implies that the 2-spinon dynamic structure factor 
vanishes in a square-root cusp: 



S$(q, 



8C ^MAm) 
tt ujfj(q) - u 2 {q) 



vM?) 



■ ui. 



(4.2) 



M^ a \q, to) vanishes more slowly, ~ [u>u(q] — lj} 1 / 2 , im- 
plying that S z a z \q,uj) drops to zero abruptly. 

At uiL(q) we find a square-root divergence (for q ^ if) 
in both the exact and the approximate transition-rate 
functions, but in the former this power-law singularity is 
accompanied by a logarithmic correction: 



/ c4(g) -"£(<?) 



111 - 



yJuj-UJ L {q)\ LO-uj L {q) 



Since the 2-spinon density of states is a step function near 

(2) 

0Jh{q), only the prefactor changes in S zz (q,uj): 



S?}(q,u) 



M(q,u) 



(4.4) 



For q — > tt the singularity at oJhiq) turns into a much 
stronger infrared singularity: 



M{tt,u)^ V2tt~C-J\yl-, (4.5) 
a; V UI 

SfX^^-fl. (4.6) 



V. SPIN AUTOCORRELATION FUNCTION 

A quantity of some interest in various experimental 
and theoretical contexts is the frequency-dependent spin 
autocorrelation function 



dte-^(S!(t)Sf). 



(5.1) 



The 2-spinon contribution to ^ zz (co) 7 

[dqSg(q, U ), 



(5.2) 



is a piecewise smooth function over the range of 2-spinon 
energies < to < tt and has singularities at u> = 0, 7r/2, tt. 
The approxi mate 2-spinon autocorrelation function in- 
ferred from ( |2.4| ) can be evaluated in terms of elliptic 
integrals. It has a step discontinuity at u> = 0, 



1 



oh, 



a logarithmic divergence at to = tt/2, 



>ir/2 



(X In ■ 



1 



\tt/2- uj\- 
and a square-root cusp at to = tt, 



(5.3) 



(5.4) 



(5.5) 



The exact 2-spinon expression has logarithmic diver- 
gences at u) — 0, tt/2, and a linear cusp at u) = tt: 



uj^O , 1 
— > oc m — . 

UJ 



(4,, *<? M »^(„ 1R _L_) 



3/2 



S^H^a (tt-c) 



(5.6) 



(5.7) 



(5.8) 



The functions and $i"^(a;) are plotted in Fig. 2. 
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VI. SUM RULES 

How important is the 2-spinon contribution to 
S zz (q 7 uj) in relation to that of other excited states? The 
key to the answer is provided by sum rules, such as the 
first frequency mo men t ( |2.2| ), which is known for all q, or 
the susceptibility (|2.3| ) , which is known for q = only, or 
the integrated intensity (static structure factor), 



oo 



(6.1) 



of which we know the grand total 
1 



It = ~ I dq I(q) = ((SI 

IT 





z\2\ 



(6.2) 



The exact 2-spinon contribution to the n th frequency mo- 
ment of S zz (q, u>), 



K n (q) 



2^ 



uJ n S zz (q,uj), 



(6.3) 



as obtained from (1.7) with ( 1 . I ) and (3.1) can be brought 
into the form 



2C r 



KW(q) = — [Ml)] n+1 kn(q), 

7T 



where 
k n (q) 



tlx 



irsinha; / , ,5 2 

2 — ( 1 — sm — tanh x 

cosh x V 2 



(6.4) 



-s(4x/tt) 



(6.5) 



For n — 2 m +1 — 1,3, ... this expression reduces to a 
polynomial in cos q, 



K (2) (a) - 

r-i / 2 \ m m 



2=0 



(-1)' 



2' 



Ki(l-cosg) m+1+( , (6.6a) 



\2l+l 



dx x(ttmhxY e . s{ixM _ 
cosh a; 



(6.6b) 



The exact sum rules for K^m+i (fl) .were shown to have 
precisely this general structure J23 : lj which, incidentally, 
is also reproduced by the frequency moments ifjim+iG?) 
of Siz (q, uj). However, the exact coefficients of the poly- 
nomial are only known for m = 0. Comparison of 



K[ 2) (q) = —k (1- cosq), k q = 0.9163..., (6.7) 

7T 



c 



with (2.2) provides one way of measuring the relative 
spectral weight of the 2-spinon excitations: 



K?\q) 
Kl(q) 



= 0.7130. 



(6.8) 



A somewhat larger share of spectral weight, K[ a \q)/ 

#1(3) = 0.8462..., is accounted for by Sj&\q,u)). 

A different way of measuring the relative 2-spinon spec- 
tral weight is provided by the static structure factor (|6.l|). 
Here, the missing spectral weight of higher-lying excita- 
tions is weighted less heavily. 

The exact 2-spinon static structure factor I^(q) = 



( 2) 

'(q) taken from (6.4) and integrated over q yields the 



total 2-spinon intensity 



4 2) = ^ / </•<■ 



cosh x 



-s(4x/tt) 



0.7289/t- (6.9) 



The total intensity of S ( z a z \q,uj) is! ~ 0.7424J T . 

The observation that si a z (q,uj) overestimates the to- 
tal 2-spinon intensity by a smaller fraction, I T °^/I^p ~ 
1.0185, than the first frequency moment of the 2-spinon 
spectral weight, K[ a \q) / K[ 2 \q) ~ 1.1868, is consistent 
with the observation made previously that it predicts 
too much spectral weight near uju(q) and too little near 
wx(«). 

At small q, where the 2-spinon continuum is very nar- 
row, all frequency moments of S^ 2 ) (q, u>) and S z a z (q, uj) 
have exactly the same ratio 



K { n a) (q) g^o 4C 
K [2 \q) ~^ * 



k = 0.8426... 



(6.10) 



(2) 

The implications of the frequency moments Kq (q) and 
(2) 

K_((q) for the singularities of the static structure fac- 
tor and the static susceptibility, respectively, will be dis- 
cussed later. 



VII. FINITE-CHAIN RESULTS 

To what extent and accuracy can the spectral - we ight 
distribution of S zz (q,u>) be reconstructed from (1.6) on 
the basis of finite-chain data for excitation energies u\ 
and transition rates M\l In a generic situation, the 
chances for succe ss m ay be remote. Convergence of the 
finite- A data for (1.6) toward the infinite- A spectral den- 
sity may only exist in an average sense, such as can be 
realized , at least in principle, by a histogram representa- 
tion of (1.6), but hardly in practice given the very coarse- 
grained spectral-weight distribution even in the largest 
systems that can be handled computationally. 

Among the ever growing collection of Bethe-ansatz 
solvable models, there exist numerous situations where 



4 



the spectral density of interest is dominated by a spe- 
cific class of excitations that can be identified in terms of 
Bethe quantum numbers. When the dynamically domi- 
nant class of excitations consists of a two-parameter con- 
tinuum, as is frequently the case, the task of reconstruct- 
ing that spectral density from finite-iV data with rea- 
sonable accuracy may be perfectly within the reach of 
state-of-the-art computational applications. 

In the case at hand, the 2-spinon excitation energies u\ 
can be evaluated for finite chains over a wide range of N 
and then again for infinite N, all via Bethe ansatz. The 
finite-iV transition rates M\ can be evaluated directly 
from the Bethe-ansatz wave function for the ground state 
and the 2-spinon states up to N = 16 and indirectly from 
the finite^Y ground-state wave function via the recursion 
methodO up to N = 28. 

The crucial point for the reconstruction of the 2-spinon 
part of the dynamic structure factor S zz (5, u>) is that 
it factorizes into two smooth functions: the density of 
states D(q,u>), which can be determined exactly via 
Bethe ansatz, and the transition rate function M(g,u>), 
toward which the finite- AT transition rates seem to con- 
verge in the following sense: pick any sequence of finite- N 
2-spinon states with energies uj\(N) and wave numbers 
qx(N) converging toward (q,u) as N — > 00. Then the 
associated scaled transition rates NM\(N) converge to- 
ward the exact transition rate function M(q,u>). 

In the main plot of Fig. 3 we show the transition rate 
functions M(tt,uj) (exact, solid line) and M^(tt, u>) (ap- 
proximate, dashed line) along with scaled finite-iV tran- 
sition rates NM\ for N = 6, 8, 28. The downward de- 
viation of AI^(n,uj) from M(tt,lo) at low frequencies is 
due to the lacking logarithmic corrections in the infrared 
divergence and the upward deviation at high frequencies 
due to the different cusp singularity at Uu(ir). 

All finite-Y data points fall close to the solid line. 
Their deviations from that line have an irregular appear- 
ance at first sight. This is attributable to the fact that an 
increasing number of spectral contributions from systems 
with increasing Y are distributed over a fixed frequency 
interval. However, when we focus on the lowest-lying ex- 
citation, for example, we see that the data points move 
away from the dashed line toward the solid line. The 
uniform convergence of this particular sequence of data 
points is best observable in the representation of the inset 
on the left of Fig. 3. 

The region near lou{^) is shown magnified in the in- 
set on the right. Here the finite- AT data converge in a 
much more complicated pattern. Nevertheless, the trend 
is clearly toward the linear behavior of the solid line and 
away from the square-root behavior of the dashed line. 

The corresponding results for q = 7r/2 are depicted in 
Fig. 4. Here the highest 2-spinon excitation for Y = 28, 
which we were unable to compute with sufficient accuracy 
via the recursion method, is not included. Even with 
the few finite-chain data points available in this case, 
the finite-size scaling behavior of the transition rates M\ 
and their convergence toward the exact transition-rate 



function is again convincingly determined. 

Given the exact asymptoti&ijinite-size gap of the lowest 
2-spinon excitation at q — irf3 



N- 



N' 



TV 

Y 



(7.1) 



1 1 

and the exact infrared divergence (4.6) of S zz (q,uj), it 
is possible to determine, under standard scaling assump- 
tions, the leading Y-dependence of the integrated inten- 
sity at q — 7T, 



/(tt.JV) 12° (In AT) 3 / 12 



(7.2) 



with to — sJlCjix. The exact coefficient, to /27t = 
0.1079..., is significantly higher than the value 0.09052 
predicted in a recent DMRG study.c3 The leading singu- 
larity of the integrated intensity for N = 00, q — ► 7r is 
then predicted to be of the form 



, . q-fir to 

I(q,ao) — ► — 

Z7T 



ln(l-i 

7T 



3/2 



(7.3) 



which is consistend with the exactly known leading 
asymptotic term of the static spin correlation functioned 
(SfSf +n ) ~ (-lJ-n-^njVa/n. 

The corresponding leading terms for the static suscep- 
tibility read: 



x^ao^^aVat, 



V~ Mtt - q) 



(7.4) 
(7.5) 
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FIG. 1. (a) Exact and (b) approximate 2-spinon dynamic 
structure factor. Both expressions are nonzero only in the 
shaded region of the (q, oj)-plane bounded by uol{o) and 
uJu(q)- 



FIG. 2. Two-spinon part of the frequency-dependent spin 
autocorrelation function. The solid line represents the ex- 
act result $i 2 ^ (u>) and the dashed line the approximate result 



FIG. 3. Two-spinon transition-rate function at q = ir. 
The solid line represents the exact result M(q,ui) and the 
dashed line the approximate result M^ a \q,uj). Also shown 
are scaled finite-chain transition rates NM\ for all 2-spinon 
excitations at q = tv of systems with N = 6, 8, 16, 28 
spins, and for lowest 2-spinon excitations also of systems with 
N = 18, 20, 26. The low-frequency and high-frequency 
parts are shown again in the insets with transformed scales 
on both axes. 



FIG. 4. Two-spinon transition-rate function at q — n/2. 
The solid line represents the exact result M(q,ui) and the 
dashed line the approximate result M' a '(g,oj). Also shown 
are scaled finite-chain transition rates NM\ for all 2-spinon 
excitations at q — n of systems with N = 8, 12, 16, 28 spins. 
The low-frequency and high-frequency parts are shown again 
in the insets with transformed scales on both axes. 
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The exact expression derived by Bougourzi, Couture, and 
Kacir for the 2-spinon contribution to the dynamic spin struc- 
ture factor S zz (q,ui) of the one-dimensional 5=1/2 Heisen- 
berg antiferromagnet at T = is evaluated for direct com- 
parison with finite-chain transition rates (TV < 28) and an 
approximate analytical result previously inferred from finite- 
TV data, sum rules, and Bethe-ansatz calculations. The 2- 
spinon excitations account for 72.89% of the total intensity 
in S zz (q,ui). The singularity structure of the exact result 
is determined analytically and its spectral-weight distribu- 
tion evaluated numerically over the entire range of the 2- 
spinon continuum. The leading singularities of the frequency- 
dependent spin autocorrelation function, static spin structure 
factor, and g-dependent susceptibility are determined via sum 
rules. 



I. INTRODUCTION 

Notwithstanding the fact that Bethe 1 found the key 
that solves the one-dimensional (ID) 5=1/2 Heisenberg 
model, 

N 

H = J ^ S; • S; + i, (1.1) 
(=1 

as early as 1931, the emergence of explicit results for 
various physical quantities (ground-state energy, 2 exci- 
tation spectrum, 3 magnetization curve, susceptibility, 4,5 
thermodynamics 6 ) was slow at first and then faster since 
around 1960. Interest in this model began to spread far 
and wide when the first compounds with quasi-lD mag- 
netic properties were synthesized and investigated exper- 
imentally. 

However, the dynamics of the ID Heisenberg antiferro- 
magnet, i.e. Hamiltonian (1.1) with J > 0, has remained 
elusive to any rigorous approach during all those years. 
An exact result for the dynamic spin structure factor 

+ oo 

S zz (q,uj) = ^ $>•"*" J dte^(Sf(t)Sf +n ), (1.2) 

-oo 

in particular, would have been of great value for the in- 
terpretation of a host of experimental data. 7 

Significant progress in the understanding of the T = 
dynamics resulted from the observation 8 that almost all 



the spectral weight in S zz (q,u>) is carried by a special 
class of Bethe-ansatz solutions with excitation energies 
(in units of J henceforth) 

u m (q) = 7rsm - cos y- —J (1.3) 

for TV —> oo, < q < ir, < q m < q. They form a two- 
parameter continuum in the (q, w)-plane bounded by the 
branches 

7r q 
w i(s)=2 sins ' wc/(?) = Trsin -. (1.4) 

These excitations were later named 2-spinon states. 
Their density of states (after rescaling by a factor 2tt/N) 
is: 8 

D( g ,q,) = e ("-^)QM«)-"). (1.5) 

The T = dynamic spin structure factor for a finite 
system with even TV and periodic boundary conditions 
can be written in the form 

S zz (q,u) = 2w^2M x 6(u-u x ), (1.6) 

A 

where M\ = \(G\S z q \X)\ 2 with S z q = N -1 ! 2 J2i e iql Sf are 
the transition rates between the singlet (St = 0) ground 
state \G) and the triplet (St = 1) states |A) with finite- 
TV excitation energies u>\. Among them are the TV (TV + 
2)/8 2-spinon excitations, which contribute most of the 
spectral weight. 

The finite-chain analysis of Ref. 8 suggested that the 
scaled transition rates TV M\ vary smoothly with q and u>. 
The consequence could be that the exact 2-spinon part of 
the dynamic structure factor is expressible, in the limit 
TV — ► oo, as a product 

S ( V(q,L0) = M(q,L0)D(q,L0), (1.7) 

with a smooth transition-rate function M(q,u>), toward 
which the scaled finite- TV transition rates converge. This 
scenario is indeed realized in the related XX model, 9,10 
where the 2-spinon density of states is given by (1.5) 
with modified spectral boundaries, and the transition- 
rate function is a constant. 8 
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II. APPROXIMATE TRANSITION RATES 

In the Heisenberg model (1.1), the finite- TV data for the 
2-spinon matrix elements indicate that M(q,u>) diverges 
at ijj = u>l(q) an d vanishes at u> = u>u(q)- In Ref. 8 the 
expression 



M {a \q,uj) 



1(3) 



(2-1) 



for the 2-spinon transition-rate function was proposed 
on the basis of this observation and the following three 
requirements: The resulting (approximate) 2-spinon dy- 
namic structure factor Si"\q,u>) = M^(q,ui)D(q,ui) 
must produce, for q = 7r, the correct infrared 
exponent. 11,12 Furthermore, it must produce the correct 
(/-dependence of the known first frequency moment, 8,13 



KM 



2tt 



uiS zz (q,uj) 



2Eg 

' 3N 



(1-cosg), (2.2) 



o 



where Eg = — TV(ln2 — 1/4) is the ground-state energy, 2 
and, via the sum rule 



x(q) = - [ —Szz(q,u), 

IT J U) 





(2.3) 



the correct value for the direct susceptibility: 4,5 x(0) 
1/-7T 2 . The resulting approximate expression, 14 



6(w - u) L (q))Q(u)u(q) - „) 



(2.4) 



for the 2-spinon dynamic structure factor has been used 
quite frequently for the interpretation of inelastic neutron 
scattering measurements on a number of quasi-lD anti- 
ferromagnets at low temperature, 7 and for comparisons 
with the results of various computational studies. 15-1 7 

It is interesting to note in this context that the ex- 
act dynamic structure factor S zz (q,u>) of the Haldane- 
Shastry model has a structure very similar to (2. 4). 18 In 
that model, as in the XX model, all the spectral weight 
of S zz (q, u>) is carried by the 2-spinon excitations. 



III. EXACT TRANSITION RATES 

A detailed assessment of the merits and limitations of 
the result (2.4) has become possible only recently through 
a remarkable new development. By approaches based on 
the concept of infinite-dimensional symmetries which had 
been developed in the context of string theory, conformal 
field theory, and quantum groups 19 Bougourzi, Couture, 



id Ka 



„20 



were able to derive the exact expression for 



the 2-spinon transition-rate function in the form 



21 



M(q,io)=-e- 



-I(t) 



where t = 2(/?i — /?2)/"" an d 
I(t)~- 



, cosh(2;c) cos(xi) — 1 „ 

dx , , , ' , — '- e x , 



x sinh(2a;) cosh x 



+ 



2 cosh j3\ 2 cosh /?2 
q = — cot _1 (sinh f3\) — cot _1 (sinh fa). 



(3.1) 

(3.2) 

(3.3a) 
(3.3b) 



By solving Eqs. (3.3) we can express the auxiliary vari- 
able t as a function of the two physical variables q,u>: 



4 , _! Ll(q)-Loi(q) 



t = —cosh 

7T 



w 2 - uj 2 L (q) 



(3.4) 



For the numerical evaluation of (3.1) we separate the sin- 
gular part from the integral (3.2): 



7(f) 

where 

h(i) = Gi(i) + h(i)-h(i) 

oo 

dx cos(xt) 



Io - In I t sinh 2 ^ ) + h(t) 



hit) 



x cosh 2 x ' 



hit) 



dx cos(xt) 
x coth 2 x ' 



l o 
Io = 7 + A(0) - / 2 (0) = 0.3677103.... 

A series expansion of s(t) = [ln^ + C — h(t)]/2, 

CO 

sm 2 (xt/2) 



(3.5) 

(3.6a) 
(3.6b) 
(3.6c) 



s(t) = I dx- 

x cosh x 



^(-l) m mln 1 + 



4m 2 



C=e Io /2 = 0.72221..., (3.7) 
brings (3.1) with t from (3.4) into closed form: 



7rf 

M(q,Lo) = Cf sinh ( — 



n 



{l + [f/(4m-2)] 2 } 
{1 + [f/4m] 2 } 2 ™ 



2\2m-l 



(3.8) 



The exact 2-spinon part of S zz (q, u>), i.e. the function 
(1.7) with the density of states (1.5) and the transition- 
rate function (3.1) evaluated numerically via (3.5) with 
(3.4) is plotted in Fig. 1(a). For comparison, the approx- 
imate result (2.4) is shown in Fig. 1(b). The two results 
look very similar, yet there are subtle differences, which 
may not matter for most experimental comparisons but 



2 



are important for comparisons with other theoretical re- 
sults. 

Both expressions diverge at the lower spectral bound- 
ary u>L(q). At the upper boundary u>u{q), Szz\q,u>) has 

(2) 

discontinuity, whereas S zz (q,ui) approaches zero contin- 
uously over a rounded shoulder. 

The structure of the exact transition rate function (3.1) 
lends itself naturally to be factorized into the approx- 
imate function (2.1) and a correction which accounts 
for the modified singularities at the boundaries of the 
2-spinon continuum: 

M(q,u>) = M {a \q,uj)^Ct/2e h ^l 2 . (3.9) 



IV. SINGULARITIES AT u L (q) AND Lou{q) 

What is the precise nature of the leading singular- 
ity in the transition-rate function M{q,u>) and in the 

(2) 

2-spinon dynamic structure factor S zz (q,ui) at the spec- 
tral boundaries u>u(q) an d u>l(q), and how do these sin- 
gularities compare with those of the approximate results 
Si a z \q, ui) and M^(q, ui)? The answer is obtained by in- 
serting (3.5) into (3.1), evaluating the leading term for 
t — ► and t —> oo, respectively, and inserting (3.4) ex- 
panded accordingly. 

At uju(q) the transition-rate function is thus found to 
approach zero linearly, 

M(q,Lo) — > ' d ^u(q) -u], (4.1) 

7T wfj{q) -w 2 L (q) 

which implies that the 2-spinon dynamic structure factor 
vanishes in a square-root cusp: 



M^ a \q,Lo) vanishes more slowly, ~ [uijj(q) — to] 1 ! 2 , im- 
plying that Si"\q,u>) drops to zero abruptly. 

At u>L(q) we find a square-root divergence (for q ^ ir) 
in both the exact and the approximate transition-rate 
functions, but in the former this power-law singularity is 
accompanied by a logarithmic correction: 



M(q,io)' 



fCj2 K(g)-o>2(g) 

7T y u> L (q) 



x ; 1 = ,/ln ] -— . (4.3) 

y/ui - u> L (q) \/ LO-LO L (q) 
Since the 2-spinon density of states is a step function near 

(2) 

uiL^q), only the prefactor changes in S zz (q,ui): 



M(q,io) 



^ 2 (g)-^(g)' 



(4.4) 



For q — ► 7r the singularity at wl(?) turns into a much 
stronger infrared singularity: 



M(tt, uj f-^ V^rC- Jin - 



V 7T Ul V Ul 



(4.5) 
(4.6) 



V. SPIN AUTOCORRELATION FUNCTION 

A quantity of some interest in various experimental 
and theoretical contexts is the frequency-dependent spin 
autocorrelation function 



+ oo 



^> zz (ui) = / dte- tut {Sf(t)Sf) 



(5.1) 



The 2-spinon contribution to <& zz (uj), 

IT 

^ 2 }n = - [dqS£Xq,u), (5.2) 



is a piecewise smooth function over the range of 2-spinon 
energies < ui < 7r and has singularities at ui = 0, tt/2, 7r. 
The approximate 2-spinon autocorrelation function in- 
ferred from (2.4) can be evaluated in terms of elliptic 
integrals. It has a step discontinuity at ui = 0, 



$W(u,)^°I + 0(u,), 

7T 

a logarithmic divergence at ui = 7r/2, 
^H^ocln 1 



|tt/2- ui\ : 
and a square-root cusp at ui = 7r, 

<J>(")(w) W -^OC ^TT - Ul. 



(5.3) 



(5.4) 



(5.5) 



The exact 2-spinon expression has logarithmic diver- 
gences at ui = 0, 7r/2, and a linear cusp at ui = 7r: 



In 



tt/2- 



3/2 



(5.6) 



(5.7) 



<Z>( 2 )(ui)^«:(ir-ui) (5.8) 
The functions <b z 2 )(u>) and <b z a )(uj) are plotted in Fig. 2. 



3 



VI. SUM RULES 

How important is the 2-spinon contribution to 
S zz {q,u>) in relation to that of other excited states? The 
key to the answer is provided by sum rules, such as the 
first frequency moment (2.2), which is known for all q, or 
the susceptibility (2.3), which is known for q = only, or 
the integrated intensity (static structure factor), 



1(2) 



2tt 



Szz(q,u), 



(6.1) 



of which we know the grand total: 



h = - 

IT 



I(q) = ((SfY) 



(6.2) 



The exact 2-spinon contribution to the n th frequency mo- 
ment of S zz (q,u>), 



K„(q) 



2tt 



u) n S zz (q,u)), 



as obtained from (1.7) with (1.5) and (3.1) can be brought 
into the form 



2C 



where 



kn(q) = dx 



Ki 2 \q) = -^Wu{q)Y + 1 k n {q), 



x sinh x j . 9 q 2 
~2 — I 1 — sin - tanh ; 
cosh x \ 2 



(6.4) 



-s(4x/tt) 



(6.5) 



For n = 2m +1 = 1,3, ... this expression reduces to a 
polynomial in cosg, 

4%M = 



7T V 2 



E 

;=o 



; K/ (l - cos g ) m+1 +', (6.6a) 



cosh x 



(6.6b) 



The exact sum rules for A'2m+i(?) were shown to have 
precisely this general structure, 22,23 which, incidentally, 
is also reproduced by the frequency moments K^+iil) 
of Si a z \q, u>). However, the exact coefficients of the poly- 
nomial are only known for m = 0. Comparison of 



K {2 \q) = -k (1 - cosg), k = 0.9163..., (6.7) 

7T 



c 



with (2.2) provides one way of measuring the relative 
spectral weight of the 2-spinon excitations: 



Ki(q) 



0.7130. 



(6.8) 



A somewhat larger share of spectral weight, K[ a \q)/ 

A' 1 (g) = 0.8462..., is accounted for by S z "\q,u>). 

A different way of measuring the relative 2-spinon spec- 
tral weight is provided by the static structure factor (6.1). 
Here, the missing spectral weight of higher-lying excita- 
tions is weighted less heavily. 

The exact 2-spinon static structure factor I^ 2 \q) = 

(2) 

Kq (q) taken from (6.4) and integrated over q yields the 
total 2-spinon intensity 



r(2) 



ac r r 2 

/ dx — e-'^l^ ~ 0.7289/t- (6.9) 
7r d J cosh x 





The total intensity of S zz \q,u) is 8 I [ T a> ~ 0.74247 T . 

The observation that S z "\q,u>) overestimates the to- 

(«)/r( 2 ) ^ 



(a) 



(6.3) tal 2-spinon intensity by a smaller fraction, Ijf'/ll 



1.0185, than the first frequency moment of the 2-spinon 
spectral weig ht, K{ a) (q)/K{ 2) (q) ~ 1.1868, is consistent 
with the observation made previously that it predicts 
too much spectral weight near u>u(q) an d too little near 
"L(q)- 

At small q, where the 2-spinon continuum is very nar- 
row, all frequency moments of S z 2 z (q,<^) and S z "\q,u>) 
have exactly the same ratio 



K ( n a) (g) 4C 
K (2 \q) ~^ * 



K = 0.8426. 



(6.10) 



(2) 

The implications of the frequency moments Kq ' \q) and 

K-I(q) for the singularities of the static structure fac- 
tor and the static susceptibility, respectively, will be dis- 
cussed later. 



VII. FINITE-CHAIN RESULTS 

To what extent and accuracy can the spectral-weight 
distribution of S zz {q,u>) be reconstructed from (1.6) on 
the basis of finite-chain data for excitation energies u>\ 
and transition rates M\? In a generic situation, the 
chances for success may be remote. Convergence of the 
finite- TV data for (1.6) toward the infinite- TV spectral den- 
sity may only exist in an average sense, such as can be 
realized, at least in principle, by a histogram representa- 
tion of (1.6), but hardly in practice given the very coarse- 
grained spectral-weight distribution even in the largest 
systems that can be handled computationally. 

Among the ever growing collection of Bethe-ansatz 
solvable models, there exist numerous situations where 
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the spectral density of interest is dominated by a spe- 
cific class of excitations that can be identified in terms of 
Bethe quantum numbers. When the dynamically domi- 
nant class of excitations consists of a two-parameter con- 
tinuum, as is frequently the case, the task of reconstruct- 
ing that spectral density from finite- TV data with rea- 
sonable accuracy may be perfectly within the reach of 
state-of-the-art computational applications. 

In the case at hand, the 2-spinon excitation energies u>\ 
can be evaluated for finite chains over a wide range of TV 
and then again for infinite TV, all via Bethe ansatz. The 
finite-TV transition rates M\ can be evaluated directly 
from the Bethe-ansatz wave function for the ground state 
and the 2-spinon states up to TV = 16 and indirectly from 
the finite- TV ground-state wave function via the recursion 
method 17 up to TV = 28. 

The crucial point for the reconstruction of the 2-spinon 
part of the dynamic structure factor S zz (q,u>) is that 
it factorizes into two smooth functions: the density of 
states D{q,u>), which can be determined exactly via 
Bethe ansatz, and the transition rate function M(q,u>), 
toward which the finite- TV transition rates seem to con- 
verge in the following sense: pick any sequence of finite- TV 
2-spinon states with energies u>\(N) and wave numbers 
q\(N) converging toward (q,u>) as TV — ► oo. Then the 
associated scaled transition rates TV Ma (TV) converge to- 
ward the exact transition rate function M{q,u>). 

In the main plot of Fig. 3 we show the transition rate 
functions M(ir, ui) (exact, solid line) and M( a \ir, ui) (ap- 
proximate, dashed line) along with scaled finite- TV tran- 
sition rates TV M\ for TV = 6,8, ...,28. The downward de- 
viation of M( a \ir, ui) from M(ir, ui) at low frequencies is 
due to the lacking logarithmic corrections in the infrared 
divergence and the upward deviation at high frequencies 
due to the different cusp singularity at uju(^)- 

All finite-TV data points fall close to the solid line. 
Their deviations from that line have an irregular appear- 
ance at first sight. This is attributable to the fact that an 
increasing number of spectral contributions from systems 
with increasing TV are distributed over a fixed frequency 
interval. However, when we focus on the lowest-lying ex- 
citation, for example, we see that the data points move 
away from the dashed line toward the solid line. The 
uniform convergence of this particular sequence of data 
points is best observable in the representation of the inset 
on the left of Fig. 3. 

The region near u>jj{^) is shown magnified in the in- 
set on the right. Here the finite- TV data converge in a 
much more complicated pattern. Nevertheless, the trend 
is clearly toward the linear behavior of the solid line and 
away from the square-root behavior of the dashed line. 

The corresponding results for q = tt/2 are depicted in 
Fig. 4. Here the highest 2-spinon excitation for TV = 28, 
which we were unable to compute with sufficient accuracy 
via the recursion method, is not included. Even with 
the few finite-chain data points available in this case, 
the finite-size scaling behavior of the transition rates M\ 
and their convergence toward the exact transition-rate 



function is again convincingly determined. 

Given the exact asymptotic finite-size gap of the lowest 
2-spinon excitation at q = 7r, 24 



n^co a 
TV' 



7T 

T 



(7.1) 



and the exact infrared divergence (4.6) of S z V(q,u), it 
is possible to determine, under standard scaling assump- 
tions, the leading TV-dependence of the integrated inten- 
sity at q = 7r, 



7( T) JV) ^ ^(lnTV) 3 / 2 



(7.2) 



with mo = y /, 2G'/ir. The exact coefficient, mo/27r = 
0.1079..., is significantly higher than the value 0.09052 
predicted in a recent DMRG study. 25 The leading singu- 
larity of the integrated intensity for TV = oo, q — ► 7r is 
then predicted to be of the form 



I(q, oo) 



q^TT mo 

' 2^ 



In 1 



IT J 



3/2 



(7.3) 



which is consistend with the exactly known leading 
asymptotic term of the static spin correlation function 26 
(SfSf +n )~(-l)»n-\hiny/yn. 

The corresponding leading terms for the static suscep- 
tibility read: 



x(tt,TV) — > — TVVTV, 



oc 



\/~ m (7r - q) 



(7.4) 
(7.5) 
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FIG. 1. (a) Exact and (b) approximate 2-spinon dynamic 
structure factor. Both expressions are nonzero only in the 
shaded region of the (q, u;)-plane bounded by iOL(q) and 

FIG. 2. Two-spinon part of the frequency-dependent spin 
autocorrelation function. The solid line represents the ex- 
act result <S>^)(uj) and the dashed line the approximate result 



FIG. 3. Two-spinon transition-rate function at q = w. 
The solid line represents the exact result M(q, to) and the 
dashed line the approximate result M^ a \q, lo). Also shown 
are scaled finite-chain transition rates NM\ for all 2-spinon 
excitations at q = w of systems with N = 6,8, 16, 28 
spins, and for lowest 2-spinon excitations also of systems with 
N = 18, 20, 26. The low-frequency and high-frequency 
parts are shown again in the insets with transformed scales 
on both axes. 

FIG. 4. Two-spinon transition-rate function at q = w/2. 
The solid line represents the exact result M(q, lo) and the 
dashed line the approximate result M^ a \q,Lo). Also shown 
are scaled finite-chain transition rates NM\ for all 2-spinon 
excitations at q = w of systems with N = 8, 12, 16, 28 spins. 
The low-frequency and high-frequency parts are shown again 
in the insets with transformed scales on both axes. 
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